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Abstract: Over almost all of minimal supergravity (mSUGRA or CMSSM) model pa- 
rameter space, there is a large overabundance of neutralino cold dark matter (CDM). We 
find that the allowed regions of mSUGRA parameter space which match the measured 
abundance of CDM in the universe are highly fine-tuned. If instead we invoke the Peccei- 
Quinn- Weinberg- Wilczek solution to the strong CP problem, then the SUSY CDM may 
consist of an axion/axino admixture with an axino mass of order the MeV scale, and where 
mixed axion/axino or mainly axion CDM seems preferred. In this case, fine-tuning of the 
relic density is typically much lower, showing that axion/axino CDM (aaCDM) is to be 
preferred in the paradigm model for SUSY phenomenology. For mSUGRA with aaCDM, 
quite different regions of parameter space are now DM-favored as compared to the case of 
neutralino DM. Thus, rather different SUSY signatures are expected at the LHC in the 
case of mSUGRA with oaCDM, as compared to mSUGRA with neutralino CDM. 
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1. Introduction 

A wide array of astrophysical data point to us living in a universe comprised of 4% baryons, 
~ 25% cold dark matter (CDM) and ~ 70% dark energy. In fact, the cosmic abundance of 
CDM has been recently measured to high precision by the WMAP collaboration]!]], which 
finds 

Q CDM h 2 = 0.110 ±0.006, (1.1) 

where £1 = p/p c is the dark matter density relative to the closure density, and h is the scaled 
Hubble constant. No particle present in the Standard Model (SM) of particle physics has 
the correct properties to constitue the CDM, so some form of new physics is needed. It is 
compelling, however, that candidate CDM particles do emerge naturally from two theories 
which provide solutions to longstanding problems in particle physics. 

The first problem- known as the gauge hierarchy problem- arises due to quadratic 
divergences in the scalar sector of the SM. These divergences lead to scalar masses blowing 
up to the highest scale in the theory (e.g. in grand unified theories (GUTS), the GUT 
scale M GU t ~ 2 x 10 16 GeV), unless an enormous fine-tuning of parameters is invoked. 
One solution to the gauge hierarchy problem occurs by introducing supersymmetry (SUSY) 
into the theory. The inclusion of softly broken SUSY leads to a cancellation of quadratic 
divergences between fermion and boson loops, so that only log divergences remain. The log 
divergence is soft enough that vastly different scales remain stable within a single effective 
theory. In SUSY theories, the lightest neutralino emerges as an excellent WIMP CDM 
candidate. Gravity-mediated SUSY breaking models (supergravity, or SUGRA) contain 
gravitinos with weak-scale masses. SUGRA models experience tension due to possible 
overproduction of gravitinos in the early universe, leading to an overabundance of CDM. 
In addition, gravitinos usually decay during or after Big Bang nucleosynthesis (BBN), 
and their energetic decay products may disrupt the successful calculations of light element 
abundances, which otherwise maintain good agreement with observation. This tension in 
SUGRA models is known as the gravitino problem. 

The second problem is the strong CP problem Q. An elegant solution to the strong 
CP problem was proposed by Peccei and Quinn (PQ) many years ago||]. The PQ solu- 
tion automatically predicts the existence of a new particle (WW)0: the axion a. While 
the original PQWW axion was soon ruled out, models of a nearly "invisible axion" were 
developed in which the PQ symmetry breaking scale was moved up to energies of order 
f a ~ 10 9 — 10 12 GeV]|, g. The axion also turns out to be an excellent candidate particle 
for CDM in the universe 0. 

Of course, it is highly desirable to simultaneously account for both the strong CP 
problem and the gauge hierarchy problem. In this case, it is useful to invoke supersymmet- 
ric models which include the PQWW solution to the strong CP problem |8|. In a SUSY 
context, the axion field is just one element of an axion supermultiplet. The axion super- 
multiplet contains a complex scalar field, whose real part is the -R-parity even saxion field 
s(x), and whose imaginary part is the axion field a(x). The supermultiplet also contains 
an i?-parity odd spin-^ Majorana field, the axino d||]. The saxion, while being an i?-parity 
even field, nonethless receives a SUSY breaking mass likely of order the weak scale. The 
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axion mass is constrained by cosmology and astrophysics to lie in a favored range 10 -2 
eV~ m a ~ 1CP 5 eV. The axino mass is very model dependent [[n], 12, |l3|| , depending 
heavily on the exact form of the superpotential and the mechanism for SUSY breaking. 
In supergravity models, it may be of order the gravitino mass m 3 / 2 ~ TeV, or as low as 
m^ 2 / f a ~keV. Conditions for realizing these extremes are addressed in [12). Here, we will 
try to avoid explicit model-dependence, and adopt m a as lying within the general range 
of keV-GeV, as in numerous previous works [pd|, |i~3| , |14| , 15, pp . An axino in this mass 



range would likely serve as the lightest SUSY particle (LSP), and is also a good candidate 
particle for cold dark matter [11, 13]. 

In a previous paper [|l6|, we investigated supersymmetric models wherein the PQ so- 
lution to the strong CP problem is assumed. For definiteness, we restricted the analysis 
to examining the paradigm minimal supergravity (mSUGRA or CMSSM) model[17|. We 
were guided in our analysis by considering the possibility of including a viable mechanism 
for baryogenesis in the early universe. In order to do so, we needed to allow for re-heat 
temperatures after the inflationary epoch to reach values Tr ~ 10 6 GeV. We found that 
in order to sustain such high re-heat temperatures, as well as generating predominantly 
cold dark matter, we were pushed into mSUGRA parameter space regions that are very 
different from those allowed by the case of thermally produced neutralino dark matter. In 
addition, we found that very high values of the PQ breaking scale f a /N of order 10 11 — 10 12 
GeV were needed, leading to the mSUGRA model with mainly axion cold dark matter, but 
also with a small admixture of thermally produced axinos, and an even smaller component 
of warm axino dark matter arising from neutralino decays. The favored axino mass value 
is of order 100 keV. We note here recent work on models with dominant axion CDM ex- 
plore the possibility that axions form a cosmic Bose-Einstein condensate, which can allow 
for the solution of several problems associated with large scale structure and the cosmic 
background radiation fl8||. 

In this paper, we will examine the mSUGRA model under the assumption 1. of neu- 
tralino CDM and 2. that mixed axion/axino DM (aaDM) saturates the WMAP measured 
abundance 1 . To compare the two DM scenarios, we will evaluate a measure of fine-tuning 
in the relic abundance 

d log at 

with respect to variations in fundamental parameters a« of the model. Such a measure 
of relic abundance fine-tuning was previously calculated in Ref. [2C] in the context of just 
neutralino dark matter. Here, we will expand upon this and also consider fine-tuning of the 
relic density in the case of mixed aaDM. Our main conclusion is that the relic abundance of 
DM is much less fine-tuned in the case of mixed aaCDM, as compared to neutralino CDM. 
Thus, we find that mixed aaCDM is theoretically preferable to neutralino CDM, at least in 
the case of the mSUGRA model, and probably also in many cases of SUGRA models with 
non-universal soft SUSY breaking terms. 

We will restrict our work to cases where the lightest neutralino Z\ is either the LSP or 
the next-to-lightest SUSY particle (NLSP) with an axino LSP; the case with a stau NLSP 

lr The possibility of mixed aaCDM was suggested in the context of Yukawa- unified SUSY in Ref. Jl9j|. 
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and an axino LSP has recently been examined in Ref. [14|. Related previous work on axino 
DM in mSUGRA can be found in Ref. p]. 

The remainder of this paper is organized as follows. In Sec. ||, we calculate the neu- 
tralino relic abundance fine-tuning parameter A ^ in the mSUGRA model due to variation 



in parameters mo and We find, in good agreement with Ref. [|20|, that the WMAP 

allowed regions are all finely-tuned for low values of tan/3. For much higher tan /3 ~ 50, the 
fine-tuning is much less with respect to mo and mj/2, but nevertheless high with respect to 
tan/3. In Sec. H we review the gravitino problem, leptogenesis and the cosmological pro- 
duction of axion and axino dark matter. In Sec. [3], we calculate the fine-tuning parameter 
A a a for mixed aaCDM under the assumption of a very light axino with rria ~ 0.1 — 1 MeV. 
The fine-tuning is always quite low, for both cases of mixed axino/axion CDM and mainly 
axion CDM. In the case of mainly axino CDM, we find the scenario less well-motivated 
since for high values of Tr ~ 10 6 GeV, the value of m„ C 0.1 MeV, making the axino 
mainly warm DM instead of cold DM. In Sec. [|, we present a summary and conclusions. 



2. Fine-tuning in mSUGRA with neutralino cold dark matter 
2.1 Overview 

We adopt the mSUGRA model @ template model for examining the issue of fine- 
tuning in cases of neutralino CDM vs. aaCDM. The mSUGRA parameter space is given 
by 

m , mi/2, A), tan/3, sign(fi), (2.1) 

where mo is the unified soft SUSY breaking (SSB) scalar mass at the GUT scale, m^ 2 
is the unified gaugino mass at Mqut, is the unified trilinear SSB term at Mqut and 
tan/3 = v u /vd is the ratio of Higgs field vevs at the weak scale. The GUT scale gauge and 
Yukawa couplings, and the SSB terms are evolved using renormalization group equations 
(RGEs) from M GUT to 

m weak} at which point electroweak symmetry is broken radiatively, 
owing to the large top quark Yukawa coupling. At m wea k, the various sparticle and Higgs 
boson mass matrices are diagonalized to find the physical sparticle and Higgs boson masses. 
The magnitude, but not the sign, of the superpotential n parameter is determined by the 
EWSB minimization conditions. 

We adopt the Isasugra subprogram of Isajet to generate sparticle mass spectrapl|. 
Isasugra performs an iterative solution of the MSSM two-loop RGEs, and includes an 
RG-improved one-loop effective potential evaluation at an optimized scale, which accounts 



for leading two- loop effects [22]. Complete one- loop mass corrections for all sparticles and 
Higgs boson masses are included p3]. For the neutralino relic density, we use the IsaReD 
subprogram of Isajet [p^]. 

Our measure of fine-tuning in the neutralino relic density, A^ , is calculated by con- 
structing a grid of points in mo — my 2 space. At each point, the change in h 2 corre- 
sponding to a change in either mo or is calculated for both a positive and negative 
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parameter change, using 



O^/i 2 dai ^zj 1 " 2 



(2.2) 



where a$ = rriQ or m\/2- For each a^, the largest A ai is selected from the results for both 
the positive and negative change. To construct the overall total , the individual values 
are added in quadrature: 



We may also consider fine-tuning due to variation in Aq and tan /3. Variation in Aq yields 
tiny variations in /i 2 unless one moves close to the stop co-annihilation region (see Fig. 
|8| in Sec. |2.3| ). Variation in tan/3 gives a slight effect on the relic density unless tan/3 



becomes very large. In this case, uia decreases[25| to the extent that txia ~ 2m^, and 
neutralino annihilation rates are greatly increased due to the ^4-resonance. Then, variation 
in tan /3 mainly shifts the location of the ^-resonance in the ttlq vs. m\/2 plane. Moving on 
and off the resonance is already accounted for by varying ijiq and m\/2- Nevertheless, in 



Sec. we present results due to including A and tan f3 in the fine-tuning calculation. 
2.2 Results from variation of rriQ and m 1( / 2 

Our first results are shown in Fig. Jlj, where we show in frame a), contours of ^l^^n 2 in the 
rriQ vs. m\/2 mSUGRA plane for Aq = 0, tan/3 = 10 and fi > 0. We also take mt = 172.6 
GeV. The well-known red regions are excluded either due to a stau LSP (left-side) or lack 
of appropriate EWSB (lower and right side). The gray-shaded region is excluded by LEP2 
chargino searches (jn^ > 103.5 GeV), and the green shaded region denotes allowable 
points with ft^/i 2 < 0.11. The region below the orange dashed contour is excluded by 
LEP2 Higgs searches, which require > 114.4 GeV; here, we actually require > 111 
GeV to reflect a roughly 3 GeV error on the RGE-improved one-loop effective potential 
calculation of m^,. 



The well-known (green-shaded) hyperbolic branch/focus point (HB/FP) region[26] 
stands out on the right side, where [i becomes small and the Z\ becomes a mixed bino- 
higgsino state. On the left edge, the very slight stau co-annihilation region p7|| is barely 
visible. We also plot contours of ^^_h 2 ranging from 5 to 80. In most of the mSUGRA 
parameter space, the relic abundance is 1-3 orders of magnitude higher than the WMAP 
measured value. The valley in Q^h 2 around m^ 2 ~ 400 GeV is due to the turn-on of the 
ZiZi — > tt annihilation mode. 

In Fig. ||, we show the neutralino relic density as a 3-d plot in the mo vs. m\/2 plane, 
to gain extra perspective. The level of fine-tuning corresponds to the slope of the surface. 
We see that in most of parameter space, the slope is relatively small, i.e. the plateau is 
nearly flat. However, in this region, the relic density is far too high. In the regions where 
^Zi^ 2 ~ 0-l> then the slope is extremely steep, corresponding to large fine-tuning: a small 
variation in fundamental parameters leads to a large change in relic density. 



2 We note that Ref. [^o| consider fine-tuning versus variation in mi, and m t . We consider these as fixed 
SM parameters, much as Mz is fixed. 
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Figure 1: In the mo vs. m\i% plane of the mSUGRA model for Ao ~ 0, tan/3 = 10 and /i > 0, we 
plot a), contours of fl^^h 2 and b). regions of fine-tuning parameter Aj; . 
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Figure 2: A 3-d plot of neutralino relic density in the mo vs. m 1 / 2 plane of the mSUGRA model 
for A = 0, tan (3 = 10 and /i > 0. 



In Fig. [T] b)., we show regions of fine-tuning parameter . A value of A^ ~ 
corresponds to no fine-tuning (a flat slope in ^j^/i 2 versus variation in all parameters), while 
higher values of A^ give increased fine-tuning in the relic density. We see immediately 
from the figure that the vast majority of parameter space, where ^^_h 2 is much too large, 
is also not very fine-tuned. However, the HB/FP region, where fi — > 0, has a very high 
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Figure 3: A plot of a), £1% h 2 and b). Ag versus too for fixed values of m±/2 = 250 GeV (blue) 
and to 1/2 = 500 GeV (red), in mSUGRA with A a = 0, tan/3 = 10 and fi > 0. 



fine-tuning, with Ag ranging from 20-100! There are also regions of substantial fine-tuning 
adjacent to the LEP2 chargino mass excluded region, due to rapid changes in ^^/i 2 as 
one approaches the Z%Z\ — > h annihilation resonance [28], and also some fine-tuning at the 
turn on of Z±Zi — > ti. Finally, we see a very narrow region of fine-tuning extending along 
the stau co-annihilation region. 

To get a better grasp, we plot in Fig. ||] a slice out of parameter space at mi/ 2 = 250 
and 500 GeV, showing in a), ^^/i 2 and in b). Ag versus mo- We see the slope in fig h 2 
is very steep in the HB/FP region, leading to Ag ~ 30 (50) for lower (higher) m\/2 values. 
In contrast, in the stau co-annihilation region, where Sl^/i 2 ~ 0.11, the value of A^ ~ 3 
(12) for low (high) m\i2- The lower m\/2 value has only moderate fine-tuning since it is 
getting close to the "bulk" annihilation region^], where Z\Z\ annihilation is enhanced 
via light t-channel slepton exhange diagrams. 

To gain a better perspective on the stau co-annihilation region, in Fig. || we show a 
blown-up portrait of the low itiq region of parameter space. The "turn-around" in the green- 
shaded WMAP allowed region in frame a), is due to the impact of the bulk annihilation 
region. Most of this area lies below the mj, = 111 GeV contour, and thus gives rise to 
Higgs bosons that are too light. In frame b). is a blow-up of the fine-tuning parameter 
A^ . We see that the major portion of the stau co-annihilation region is fine-tuned, with 
the possible exception of the region lying just below the LEP2 rrih bound, where mixed 
bulk/co- annihilation occurs. 

In Fig. ||, we show contours of relic density and Ag for tan/3 = 30. At higher values 
of tan/3, the b and r Yukawa couplings increase in magnitude, and enhance neutralino 
annihilation into bb and rf final states. Overall, we see a similar picture to that shown 
in Fig. |H in that the HB/FP region has extreme fine-tuning of the relic density, while the 
stau co-annihilation region is also fine-tuned, but somewhat less so. The bulk annihilation 
region is again excluded by the LEP2 rrih bound. 

We plot in Fig. n the mSUGRA plane for tan /3 = 53. In this CclSG, Si large new green- 
shaded region is opening up along the low mo edge of parameter space. This is due to three 
effects occuring at large tan/3. 1. The f\ mass decreases with tan/3, leading to increased 
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Figure 4: A plot of a), contours of f^/i 2 anc ^ i n mSUGRA with Aq = 0, tan/? = 10 and 

/i > 0. This plot zooms in on the stau co-annihilation region. 
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Figure 5: In the tuq vs. vn\i% plane of the mSUGRA model for Aq — 0, tan f3 — 30 and /i > 0, we 
plot a), contours of ft^^h 2 and b). regions of fine-tuning parameter A^. 



annihilation into tt final states; this increases the area of the bulk annihilation region. 2. 
The tau and b Yukawa couplings f T and ft, increase, thus enhancing annihilation into rf and 
bb final states. 3. The value of txia is decreasing while the width T^ is increasing (due to 
increasing Yukawa couplings that enter the A decay modes), so that Z\Z\ — > A^*' — > bb, rf 
increases: i.e. we are entering the A resonance annihilation region [ 30 1, which enhances the 
neutralino annihilation cross section in the early universe, thus lowering the relic density. 
In the case of tan f3 = 53, we see that the HB/FP region is still highly fine-tuned. However, 
broad portions of the low m u mSUGRA parameter space around m^j ~ 300 — 600 have 
Aj: ~ 3 due to an overlap of bulk annihilation through staus, stau co-annihilation and A- 
resonance annihilation. Another low fine-tuning and relic density consistent region occurs 
at mxji ~ 1200 GeV, where one sits atop the ^-resonance. 
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Figure 6: In the mo vs. m\i% plane of the mSUGRA model for Aq — 0, tan /3 — 53 and /i > 0, we 
plot a), contours of fl^^h 2 and b). regions of fine-tuning parameter Aj; . 
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Figure 7: In the mo vs. m 1 / 2 plane of the mSUGRA model for Ao — 0, tan/3 = 55 and fi > 0, we 
plot a), contours of Clg h and b). regions of fine-tuning parameter Ag . 

In Fig. ^, we show the mSUGRA plane for tan (3 = 55. Here, the ^4-resonance annihi- 
lation region is fully displayed, and the A width is even larger. While much of the HB/FP 
region is still very fine-tuned, the regions of annihilation though the broad A resonance 
yield relatively low fine-tuning, especially if one sits right on the resonance, or sits in the 
resonance/bulk/co- annihilation overlap region at low m,Q and low m\i2- 

2.3 Results from variation in Aq and tan/3 

As mentioned earlier, including Aq into the measure of fine-tuning typically yields only 
a small effect, unless one is near the top-squark co-annihilation region. This is because 
variation in Aq mainly leads to different mixing in the third generation scalar system, and 
for most of mSUGRA parameter space, affects mainly the top squark mass eigenstates. To 
show this explicitly, we plot in Fig. ga). the value of Qj? h 2 and in frame b). the value of 
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Figure 8: Plot of a). h 2 and b). \Aa \ versus Aq for two slices out of mSUGRA parameter 
space: 1. m = 1.5 TeV and m 1 / 2 = 250 GeV (blue curves) and 2. m = 2 TeV and = 750 

GeV (red dashed curve), for tan/3 = 10 and fi > 0. 




Figure 9: Plot of a). Cl^h? and b). \ Atan p | versus tan for two slices out of mSUGRA parameter 
space: 1. mo = 1.5 TeV and = 250 GeV (blue curves) and 2. mo = 2 TeV and = 750 

GeV (red dashed curve), for Aq = and /i > 0. 



|A^ | versus variation in Aq for two cases: 1. mo = 1.5 TeV and m l j 2 = 250 GeV (blue 
curves) and 2. mo = 2 TeV and mi/2 = 750 GeV (red dashed curve), for tan/3 = 10 and 
fi > 0. In these two cases, the slope of ^■^ i h 2 is rather mild, leading to a contribution to 
|A^ | of typically 1 or less. The exception comes for the m^ = 750 GeV curve around 
Aq ~ —4 TeV, where indeed the value of m^ is rapidly becoming lighter, and feeding into 
the relic density calculation. 

In Fig. H we show a), the relic density and b). |A tan ^| versus variation in tan/3 for 
1. itiq = 1.5 TeV and ni\/2 = 250 GeV (blue curves) and 2. tuq = 2 TeV and m,^ = 750 
GeV (red dashed curve), for Aq = and fi > 0. For most of the tan/3 values, the relic 
density varies only slowly with tan /3, leading to only small contributions to A. When tan /3 
becomes of order 50, then wia is rapidly decreasing, and Ta is rapidly increasing, leading 
to a high rate of neutralino annihilation through the ^4° resonance. In this case, while 
fine-tuning with respect to tuq and m 1( / 2 is low, fine-tuning with respect to tan/3 is high. 

In Fig. 1C, we show the value of A including contributions from variation in mo, my 2 
and tan/3, for the large values of a), tan/3 = 53 and b). tan/3 = 55. Here, over essentially 
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Figure 10: Plot of A in the mo vs. m\i% plane including variation of mo, m\i% and tan/3 for 
A = 0, \i > and a), tan/3 = 53 and b). tan/3 = 55. 

all of parameter space, the value of A has increased to much larger values than those for 
A ^ as shown in Fig's [] and [7|. Thus, inclusion of tan (3 in the fine-tuning calculation shows 

that large values of tan (3 ~ 50 result in large fine-tuning of the relic density. 

3. The gravitino problem, leptogenesis, and the re-heat temperature 

In this section, we review the gravitino problem, baryogenesis via leptogenesis, and produc- 
tion of mixed axion/axino dark matter in the early universe. The reader who is familiar 
with these issues may proceed directly to Sec. ||; others may wish to follow the brief 
treatment given here and in Ref. J16| . 

3.1 The gravitino problem 

In supergravity models, supersymmetry is broken via the superHiggs mechanism. The 
common scenario is to postulate the existence of a hidden sector which is uncoupled to the 
MSSM sector except via gravity. The superpotential of the hidden sector is chosen such 
that supergravity is broken, which causes the gravitino (which serves as the gauge particle 
for the superHiggs mechanism) to develop a mass m 3 / 2 ~ rn 2 /Mpi ~ m wea k- Here, m is a 
hidden sector parameter assumed to be of order 10 11 GeV. 3 In addition to a mass for the 
gravitino, SSB masses of order m wea k are generated for all scalar, gaugino, trilinear and 
bilinear SSB terms. Here, we will assume that 7713/2 is larger than the lightest MSSM mass 
eigenstate, so that the gravitino essentially decouples from all collider phenomenology. 

In all SUGRA scenarios, a potential problem arises for weak-scale gravitinos: the 
gravitino problem. In this case, gravitinos G can be produced thermally in the early 
universe (even though the gravitinos are too weakly coupled to be in thermal equilibrium) 
at a rate which depends on the re-heat temperature Tr of the universe. The produced 
G can then decay to various sparticle-particle combinations, with a long lifetime of order 

3 In Ref. j3l]], a link is suggested between hidden sector parameters and the PQ breaking scale f a - 
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1 — 10 5 sec (due to the Planck suppressed gravitino coupling constant). The late gravitino 
decays occur during or after BBN, and their energy injection into the cosmic soup threatens 
to destroy the successful BBN predictions of the light element abundances. The precise 
constraints of BBN on the gravitino mass and Tr are presented recently in Ref. |32| . One 
way to avoid the gravitino problem in the case where m.3/2 ~ 5 TeV is to maintain a value of 
Tr ~ 10 5 GeV. Such a low value of Tr rules out many attractive baryogenesis mechanisms, 
and so here instead we assume that m 3 / 2 ~ 5 TeV. In this case, the G is so heavy that 
its lifetime is of order 1 sec or less, and the G decays near the onset of BBN. In this case, 
values of Tr as large as 10 9 GeV are allowed. 

In the simplest SUGRA models, one typically finds tuq = 7713/2- For more general 
SUGRA models, the scalar masses are in general non-degenerate and only of order m 3 / 2 p3| . 
Here for simplicity, we will assume degeneracy of scalar masses, but with thq <C m 3/ / 2 - 

3.1.1 Leptogenesis 

One possible baryogenesis mechanism that requires relatively low Tr ~ rriweak is elec- 
troweak baryogenesis. However, calculations of successful electroweak baryogenesis within 

the MSSM context seem to require sparticle mass spectra with rrih ~ 120 GeV, and 

< 



125 GeV[34|. The latter requirement is difficult (though not impossible) to achieve 
in the MSSM, and is also partially excluded by collider searches for light top squarks p5| . 
We will not consider this possibility further. 

An alternative attractive mechanism- especially in light of recent evidence for neutrino 
mass- is thermal leptogenesis . In this scenario, heavy right-handed neutrino states iVj 
{i = 1 — 3) decay asymmetrically to leptons versus anti-leptons in the early universe. The 
lepton-antilepton asymmetry is converted to a baryon-antibaryon asymmetry via sphaleron 
effects. The measured baryon abundance can be achieved provided the re-heat temperature 
Tr exceeds ~ 10 9 GeV [37]. The high Tr value needed here apparently puts this mechanism 



into conflict with the gravitino problem in SUGRA theories. 

A related leptogenesis mechanism called non-thermal leptogenesis invokes an alter- 
native to thermal production of heavy neutrinos in the early universe. In non-thermal 
leptogenesis [|3§|] , it is possible to have lower re-heat temperatures, since the iVj may be 
generated via inflaton decay. The Boltzmann equations for the B — L asymmetry have 



been solved numerically in Ref. [39|. The B — L asymmetry is then converted to a baryon 
asymmetry via sphaleron effects as usual. The baryon-to-entropy ratio is calculated in p9| , 
where it is found 



.s 



- 8.2 x 10-- x ( ) ( ^ ) ( 7 J1^ 7 ) . (3 .l 



10 6 GeV / \ niA, V0.05 eV 



where is the inflaton mass and 5 e ff is an effective CP violating phase which may be of 
order 1. Comparing calculation with data (the measured value of hr/s ~ 0.9 x 10 -10 ), a 
lower bound Tr ~ 10 6 GeV may be inferred for viable non-thermal leptogenesis via inflaton 
decay. 



A fourth mechanism for baryogenesis is Affleck-Dine poj] leptogenesis [41]. In this ap 



proach, a flat direction fa = (2H£i) 1 ^ 2 is identified in the scalar potential, which may have 
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a large field value in the early universe. When the expansion rate becomes comparable 
to the SSB terms, the field oscillates, and since the field carries lepton number, coherent 
oscillations about the potential minimum will develop a lepton number asymmetry. The 
lepton number asymmetry is then converted to a baryon number asymmetry by sphalerons 
as usual. Detailed calculations^] find that the baryon-to-entropy ratio is given by 

s 23 m v M 2 pi v 1 

where {H) is the Higgs field vev, m v is the mass of the lightest neutrino and Mpi is the 
Planck scale. To obtain the observed value of ns/s, values of Tr ~ 10 6 — 10 8 are allowed 
for m u ~ 1(T 9 - 1(T 7 eV. 

Thus, to maintain accord with either non-thermal or Affleck-Dine leptogenesis, along 
with constraints from the gravitino problem, we will aim for aaDM scenarios with Tr ~ 
10 6 - 10 8 GeV. 

3.2 Mixed axion/axino dark matter 
3.2.1 Relic axions 

Axions can be produced via various mechanisms in the early universe. Since their lifetime 
(they decay via a — > 77) turns out to be longer than the age of the universe, they can 
be a good candidate for dark matter. As we will be concerned here with re-heat temper- 
atures Tr ~ 10 9 GeV < fa/N (to avoid overproducing gravitinos in the early universe), 
the axion production mechanism relevant for us here is just one: production via vacuum 
mis-alignment . In this mechanism, the axion field a(x) can have any value ~ f a at tem- 
peratures T 3> Aqcd- As the temperature of the universe drops, the potential turns on, and 
the axion field oscillates and settles to its minimum at —Of a /N (where = 9 + arg{det m q ), 
9 is the fundamental strong CP violating Lagrangian parameter and m q is the quark mass 
matrix). The difference in axion field before and after potential turn-on corresponds to the 
vacuum mis-alignment: it produces an axion number density 

n a {t) ~ ^m a (t)(a 2 (t)), (3.3) 

where t is the time near the QCD phase transition. Relating the number density to the 
entropy density allows one to determine the axion relic density today 0: 

1 / 6 x 10-6 eV \ ^ „ 1 ( f a ,N \ ^ 
u - n ~ 4 [ ma ) ^ " 4 VlO 12 GeVy* 4 ' {6) 

where 6{ is the initial vacuum mis-alignment angle, with — 7r ~ 6i ~ n. An error estimate 
of the axion relic density from vacuum mis-alignment is plus-or-minus a factor of three. 
Axions produced via vacuum mis-alignment would constititute cold dark matter. However, 
in the event that (a 2 (t)) is inadvertently small, then much lower values of axion relic density 
could be allowed. Additional entropy production at t > tQCD can also lower the axion relic 
abundance. Taking the value of Eq. ( |3.4j ) literally, along with 0j ~ 1, and comparing 
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to the WMAP5 measured abundance of CDM in the universe, one gets an upper bound 
f a /N ~ 5 x 10 11 GeV, or a lower bound m a ~ 10~ 5 eV. If we take the axion relic density a 
factor of three lower, then the bounds change to f a /N ~ 1.2 x 10 12 GeV, and m a ~ 4 x 10~ 6 
eV. 

3.2.2 Axinos from neutralino decay 

If the a is the lightest SUSY particle, then the Z\ will no longer be stable, and can decay via 
Z\ — > 07. The relic abundance of axinos from neutralino decay (non-thermal production, 
or NTP) is given simply by 

n NTP h 2 = rn ±Q h 2 (35) 
^1 

since in this case the axinos inherit the thermally produced neutralino number density. 
The neutralino-to-axino decay offers a mechanism to shed large factors of relic density. For 
a case where mj; ~ 100 GeV and £1% h 2 ~ 10 (as can occur in the mSUGRA model at 
large itiq values) an axino mass of less than 1 GeV reduces the DM abundance to below 
WMAP-measured levels. 

The lifetime for these decays has been calculated, and it is typically in the range of 
t{Z\ — > (17) ~ 0.01 — 1 sec|l3|]. The photon energy injection from Z\ — > 07 decay into 
the cosmic soup occurs typically before BBN, thus avoiding the constraints that plague 
the case of a gravitino LSP[32|. The axino DM arising from neutralino decay is generally 



considered warm or even hot dark matter for cases with m„ ~ 1 — 10 GeV [[42]] . Thus, in the 
mSUGRA scenario considered here, where m« ~ 1 — 10 GeV, we usually get warm axino 
DM from neutralino decay. 

3.2.3 Thermal production of axinos 

Even though axinos may not be in thermal equilibrium in the early universe, they can still 
be produced thermally via scattering and decay processes in the cosmic soup. The axino 
thermally produced (TP) relic abundance has been calculated in Ref. Jl^, ft3| ], and is given 
in Ref. |^3| using hard thermal loop resummation as 

n r^ 5 ,^<^)(^) 2 (^-)(^) (3,, 

where g s is the strong coupling evaluated at Q = Tr and N is the model dependent color 
anomaly of the PQ symmetry, of order 1. For reference, we take g s (Tu = 10 6 GeV) = 0.932 
(as given by Isajet 2-loop g s evolution in mSUGRA), with g s at other values of Tr given 
by the 1-loop MSSM running value. The thermally produced axinos qualify as cold dark 
matter as long as ~ 0.1 MeV[13, 43]. 



4. Fine-tuning in mSUGRA with mixed axion/axino CDM 

In this section, we calculate the fine-tuning parameter for the dark matter relic density 
in models with mixed aaDM: A a „. Contributions to A a a are calculated from both the 
axion relic density and the thermally produced axino relic density. We do not include the 
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non-thermally produced axino relic density as it makes a tiny contribution to the total for 
the values of nij ~ 1 MeV considered here (see Figs. 2 and 3 of Ref. fl6|]). We take the 
total relic density to be 



n a ah z = n a h z + nl^h z (4.i) 

1 ( fa/N \ 7/ \ 2 6l ( 1.211 \( 10 11 GeV \ 2 / m- a \( T R \ 

(4.2) 

and calculate the total A aa exactly by differentiating (^) with respect to f a /N , Tr and 
m„. 4 We find: 

_ T R an aa h 2 _ T R e r. r fi 1n / 1.211 \ / 10 11 GcV \ 2 / m a w 1 \ I a o\ 
A - rn- a dQ. a& h 2 _ "^a e c„6 1n / l.21l \ / lO 11 GeV \ 2 / 1 \ f (4 4 \ 



and 



fa/N on a qh 2 

^a'ah 2 df a /N 



A /a/iV = o ^ ^ m7 ( 4 - 5 ) 



fa/N 
J L aa' L 



7 ( 1 



7/6 

(/ a /A0 1/6 # 2 



24 V 10 12 GeV / 

1.211 \ , n „ Tr ,2 / 1 V/ "la \/ Tfl 



ns ^) (10 n Gevf (_i_y 



0.1 GeV/ V 10 4 GeV 



(4.6) 



and 



[9i df2 a5 /i 2 _ O a /i 2 



The total fine tuning parameter is then given by 



A a , = ^A 2 fl +A 2 rta + A 2 a/iV + A 2 i . (4.8) 
We plot our first results in the f a /N vs. Tr plane, keeping m a fixed at 1 MeV: see 



Fig. 11. In frame a)., we show contours of Cl aa h 2 . The green region gives Q a a < 0.11, and 
so is consistent with WMAP. In frame b)., we show regions of fine-tuning A aa . The scale 
is shown on the right edge of the plot. Note in this case the entire plane has A aa < 2.5, so 
there is very little fine-tuning across the entire plane of parameter space. The left region, 



4 Here, one objection may be that the value of Tr does not appear as an explicit Lagrangian parameter. 
However, in the standard inflationary cosmology, the reheat temperature is related to the inflaton decay 
width via T R ~ {Z/-K 3 ) l/4 'g- 1/ \MpiT^ 1 l' 2 ^, where I> depends on the inflaton mass and couplings to 
matter. In this case, a detailed model including the inflaton field cj> would provide Tr in terms of inflaton 
Lagrangian parameters. We do not wish to bring such model-dependence into our calculations, so instead 
just adopt the value of Tr as a fundamental parameter. Also, the value of mj will appear as a Lagrangian 
parameter in the weak scale effective Lagrangian, after the effects of SUSY breaking and PQ breaking are 
taken into account. 



-14- 




S.Q S.5 10,0 10.S 11. D 11.5 12.0 12.5 

LogL/JN (GeV)] 




9.0 3.5 10.0 1D.5 11.0 11.5 12.0 12.5 

Log^/N (GeV)] 



— 12.4 
— 12.2 
2 

1.8 
1.6 
1.4 



Figure 11: Plot of a), contours of axion/axino relic density fl a 5,h 2 in the f a /N vs. Tr plane and 
b). regions of fine-tuning A a for fixed axino mass ma = 1 McV. The green region in a), has 
Qaah 2 < 0.11. The green dashed line in frame b). is where £l a d,h 2 = 0.11. 



color-coded dark blue, is the region of dominantly thermally produced axino CDM, whilst 
the right-most region, color-coded lighter blue, is dominantly axion CDM. In this region, 
the fine-tuning parameter A aa — 2.3. The intermediate region, shaded by yellow and 
purple bands, is the region of mixed aaCDM: the purple band has very low fine-tuning, 
with A a a < 1-4. The contour where mixed aaCDM saturate the WMAP measured value 
is shown by the green dashed line. The region where the highest values of Tr are found 
coincide with the region of lowest fine-tuning, with a nearly equal mix of axion and axino 
CDM. 



To gain additional perspective, in Fig. 12 we show the mixed axion/axino relic den- 
sity as a 3-d plot in the f a /N vs. Tr plane for m a = 1 MeV. The level of fine-tuning, 
corresponding to the slope of the surface, is rather low throughout, since there are no re- 
gions with a steep slope. The fine-tuning is minimal along the trough running through the 
right-center of the plot. 

To better understand the situation with mixed aaCDM, we show in Fig. ^ a slice of 
our Fig. [ll| with constant Tr = 10 5 GeV. In frame a)., we see that the value of £l a ah 2 
initially drops as f a /N increases. This is in the region of dominant axino CDM, and 
increasing f a /N decreases the axino coupling strength, and hence suppresses its thermal 
production in the early universe. As f a /N increases further, the relic abundance of axions 
steadily increases, until around f a /N ~2x 10 11 GeV there is an upswing in the relic 
abundance. This is the stable fine-tuning region since small fluctuations of parameters 
about this point do not substantially alter the axino/axino relic density. The fine-tuning 
parameter A aa is shown in frame b).. Here, we see that the fine-tuning is slightly high 
in the region of mainly axino CDM, with low f a /N, but reaches a minimum at the point 
of equal admixture. The value of A aa doesn't extend all the way to zero, in spite of the 
zero slope shown, because A aa still varies with m a , Tr and 0j. The fine-tuning parameter 
increases to the analytic value of 2.3 as f a /N increases further, into the region of mainly 
axion CDM. 
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Figure 12: A 3-d plot of axion/axino relic density in the f a /N vs. Tr plane of the ad augmented 
mSUGRA model. 
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Figure 13: Plot of a), axion/axino relic density £l a ah 2 and b). fine-tuning parameter A a a versus 
f a /N for fixed axino mass ma = 1 MeV and fixed re-heat temperature Tr = 10 5 GeV. 



A similar plot to Fig. 11 is shown in Fig. [LJ, but in this case taking m a = 0.1 MeV. 

below which 



Note that this value yields the approximate dividing line given in Refs. [13 



the thermally produced axinos would be mainly warm DM instead of cold DM. In any case, 
in frame a)., we see that the WMAP allowed region has expanded, and now values of Tr as 
high as 5 x 10 6 GeV are allowed, making the scenario consistent with at least non-thermal 
leptogenesis. The region of maximal Tr also coincides with the region of least fine-tuning, 
with a roughly equal admixture of axion and thermally produced axino DM. 

In Fig. 15 we show the contours of relic density Q aa h 2 in the m a vs. Tr plane for fixed 
value of f a /N = 4.88 x 10 11 GeV. The large value of f a /N yields a scenario with mainly 
axion CDM when the WMAP measured abundance is saturated. The green shaded region 
in frame a), is WMAP-allowed. The red dashed line shows the approximate dividing line 
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Figure 14: Plot of a), contours of axion/axino relic density fl a 5,h 2 in the f a /N vs. Tr plane and 
b). regions of fine-tuning A a a for fixed axino mass mg = O.i MeV. The green region in a), has 
Qaah 2 < 0.11. The green dashed line in frame b). is where Q a hh 2 = 0.11. 



% - 

i — i l-j 

o 







Axino 




WDM 




a) 





-6 -5 -4 -3 -2 -I 

Logfm^ino (GeV)] 




-6 -5 -4 -3 -2-10 ] 

Log[m axino (GeV)] 



2.4 
2.2 



1.8 
1.6 
1.4 



Figure 15: Plot of a), contours of axion/axino relic density fl a ah 2 in the ma vs. Tr plane and b). 
regions of fine-tuning A a a for fixed f a /N = 4.88 x 10 11 GeV (which gives mainly axion CDM along 
the line of ft a a = 0.11). The green region in a), has Sl a a,h 2 < 0.11. The green dashed line in frame 
6). is where Q a a,h 2 — 0.11. The region to the left of red-dashed line gives thermally produced warm 
axino dark matter. 



between warm and cold thermally produced axinos. In this case, the demarcation line is 
largely irrelevant, since if ^l a ah 2 — 0.11, almost all the DM is composed of cold axions, and 
a tiny admixture of warm axinos would be allowed. In frame b)., we show the regions of 
fine-tuning A a a- Since the WMAP-allowed region coincides with mainly axion CDM, the 
fine-tuning along the green dashed line is always low: A a a ~ 2.3. Note that in the scenario 
with mainly axion CDM, the value of Tr can easily reach to well over 10 7 GeV, allowing 
for non-thermal or Affleck-Dine leptogenesis. 



In Fig. 16 we show the rria vs. Tr plane for f a /N = 3 x 10 11 GeV, which gives roughly 
an equal admixture of axion and thermally produced axino DM. In this case, the value of 
Tr reaches beyond 10 8 GeV, although for very low values of nij where it is expected that 
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Figure 16: Plot of a), contours of axion/axino relic density ft a 5,h 2 in the mj vs. Tr plane and b). 
regions of fine-tuning A a j for fixed f a /N = 3 x 10 11 GeV (which gives a 50-50 mix of axion/axino 
DM along the line of ft a a — 0.11). The green region in a), has tt a 5,h 2 < 0.11. The green dashed 
line in frame b). is where ftaah 2 — 0.11. The region to the left of red-dashed line gives thermally 
produced warm axino dark matter. 



the axino will be warm DM. For this scenario, it is unclear how much mixture of warm and 
cold dark matter is cosmologically allowed. To answer the question, the velocity profile of 
the warm axinos would have to be fed into n — body simulations of large scale structure 
formation, to see how well such a mixed warm/cold DM scenario fits the data. At present, 
we are unaware of such studies. In frame &)., we see that the line of WMAP-saturated 
abundance lies nearly on top of the region of lowest fine-tuning, with A a a < 1.4. 

In Fig. 17, we show again the m„ vs. Tr plane, but this time for f a /N = 1 x 10 11 GeV, 
which gives mainly thermally produced axino DM. In this case the region to the left of the 
red-dashed line should likely be disallowed, since the dominant form of DM will be warm, 
rather than cold. The region to the right of the ma = 0.1 MeV line, in the WMAP-allowed 
region, only allows for Tr to reach a max of 10 6 GeV. Furthermore, from frame b)., we 
see that the fine-tuning parameter in this case for the WMAP-saturated region along the 
green dashed line is somewhat higher, reaching A 

ad ^ 



5. Summary and conclusions 

In this paper, we have examined the fine-tuning associated with the relic density of dark 
matter in the minimal supergravity model. We have calculated a measure of fine-tuning 
assuming two scenarios for SUSY dark matter: 1. neutralino dark matter with fine-tuning 
parameter Aj; , and 2. mixed axion/axino dark matter with fine-tuning parameter A aa . 

In the case of neutralino dark matter, we find that the WMAP-allowed regions of 
mSUGRA such as the stau co-annihilation region, the HB/FP region and the light Higgs 
/i-resonance annihilation region, are all rather highly fine-tuned, especially the HB/FP 
region, where ranges from 20-100. Only mild fine-tuning is found in the low mo, low 
m i/2 region where stau co-annihilation and bulk annihilation through t-channel slepton 
exchange overlap. If one moves to large tan (5 ~ 50, then larger regions of parameter space 
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Figure 17: Plot of a), contours of axion/axino relic density ft a a,h 2 in the ma vs. Tr plane and b). 
regions of fine-tuning A a j for fixed f a /N = 1 x 10 11 GeV (which gives mainly axino DM along the 
line of ft a a = 0.11). The green region in a), has fl a 5,h 2 < 0.11. The green dashed line in frame b). 
is where fl a d,h 2 — 0.11. The region to the left of red-dashed line gives thermally produced warm 
axino dark matter. 



which are consistent with WMAP occur. These large tan j3 regions have modest fine-tuning 
versus thq and m^, but very large fine-tuning versus tan/3. 

If instead we assume that dark matter is composed of an axion/axino admixture, rather 
than neutralinos, then we find that the relic density fine-tuning parameter is generically 
much lower: A aa ~ 1.3 — 2.5 throughout parameter space. Here, we have assumed the 
existence of a light axino with mass ra„ ~ keV-MeV. Such a light axino opens up all 
of mSUGRA parameter space to being WMAP allowed, since now the neutralino decays 
via Z\ — > a^. If the DM is dominated by thermally produced axinos, then the re-heat 
temperature Tr is generally lower than 10 6 GeV unless the axinos are actually warm dark 
matter (m a ~ 100 keV), so this scenario seems rather unlikely. However, if the PQ breaking 
scale fa/N is large, then the DM can be either a nearly equal axion/axino admixture, in 
which case fine-tuning is lowest (A aa ~ 1.3), or a dominantly axion mixture (in which case 
A a a ~ 2.3). Either scenario easily admits Tr > 10 6 GeV, which can allow for non-thermal 
leptogenesis to occur. 

The consequences of the mixed aaCDM scenario for future dark matter searches is 
as follows. For collider searches, we expect much the same collider signatures as in the 
mSUGRA model with neutralino dark matter, since we assume the Z\ is the NLSP, and 
decays far outside the collider detectors. However, all of mSUGRA parameter space is now 
WMAP-allowed, instead of just the special co-annihilation, HB/FP region and resonance 



annihilation regions. As shown in Ref. [16], the regions of WMAP-allowed neutralino 



CDM yield the lowest values of Tr, and so the stau co-annihilation, HB/FP region and h 
resonance annihilation regions are most dis-favored for the case of mixed aaCDM. 

As far as WIMP searches go, in the mixed aaCDM scenario, we expect no positive 
signals if rrig > m a . If m a > mg^, then the Z\ would still be stable (assuming i?-parity 
conservation) and WIMP direct and indirect detection signals are still possible [EBI]. In the 
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case of large axion relic abundance, which appears to us to be the favored scenario, then 
a positive signal at relic axion search experiments such as ADMX might be expected p6| , 
although solar axion searches are less likely to achieve positive results, since large values 
of fa/N are favored, leading to small axion/axino couplings. 

Our analysis has been based on the admittedly subjective basis of fine-tuning of the 
relic density of dark matter relative to model input parameters. We note here that the 
mSUGRA model already needs substantial fine-tuning in the electroweak sector in order 
to accomodate the relatively light Z boson mass in the face of limits on the soft SUSY 
breaking parameters^] (the little hierarchy problem). Our philosophy here is that less 
fine-tuning is better, and high fine-tuning in one sector is better than high fine-tuning in 
two sectors, e.g. electroweak and dark matter sectors. 

While our analysis has been restricted to the mSUGRA SUSY model, one might ask 
how general our conclusions might be. In SUSY models based on gravity mediation, with a 
neutralino LSP, the DM relic density is generically too high unless some special mechanism 
is acting to enhance the neutralino annihilation cross section in the early universe. 5 For 
instance, in SUSY models with non-universality, instead of stau or stop co-annihilation, 
one might have sbottom or sneutrino co-annihilation, or bino-wino co- annihilation: in any 
case, the mass gap between co-annihilating particles must be fine-tuned to obtain agreement 
with the measured dark matter abundance. In non-universal models with a well-tempered 



neutralino\ 48], where the neutralino bino-higgsino or bino-wino composition is adjusted to 
fit the measured relic density, other parameters (Higgs soft masses, gaugino masses) must 
be fine-tuned to get just the right "tempering", as occurs in the mSUGRA HB/FP region. 
In other models, Higgs soft mass terms can be adjusted to allow 2m % to sit atop the A 
resonance; but again, in this case, parameters must be fine-tuned (unless tan /3 is large, 
which also occurs in mSUGRA). The case where the SUSY neutralino abundance is not 
fine-tuned has long been noted: it is where squarks and sleptons are so light that t-channel 
annihilation channels are large. However, LEP2 search limits now essentially exclude all 
these regions. Thus, although we restrict our analysis here to the mSUGRA model, we 
feel this model provides a sort of microcosm for general SUSY models, in that it illustrates 
many of the features common to all SUSY models. 

Our main conclusion is this. In the world HEP community, a tremendous effort is 
underway to explore for WIMP cold dark matter, based partly on the view that SUSY 
models naturally give rise to the "WIMP-miracle" , and an excellent WIMP candidate for 
CDM. We have shown here that at least for the paradigm SUSY model- mSUGRA- usually 
a large overabundance of neutralino CDM is produced, unless one lies along a region of 
very high fine-tuning, where a slight change in model parameters leads to a large change in 
relic density: this equates to a high degree of relic density fine-tuning. Alternatively, if one 
assumes the PQWW solution to the strong CP problem within SUSY models, and a very 
light axino with rria of the order of MeV, then along with an elegant solution to the strong 
CP problem, one obtains a mixed axion/axino relic density with much less fine-tuning. 
Given our results, we would advocate that a much increased share of HEP resources be 



5 Discussion on numerous different SUGRA models with non-universality has been explored in Ref. [ H . 
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given to relic axion searches, where the global search effort has been much more limited. 



Acknowledgments 

We thank H. Summy for discussions. This research was supported in part by the U.S. 
Department of Energy. 



References 



J. Dunkley et al. [WMAP Collaboration], Astrophys. J. Suppl. 180 (2009) 306 



For a recent review, see P. Sikivie, [hep-ph/0509198| ; M. Turner, \Phys. Rept. 197 (1990) 67 
J. E. Kim and G. Carosi, arXiv:0807.3125 (2008). 



R. Peccei and H. Quinn, \Phys. Rev. Lett. 38 (1977) 1440| and \Phys. Rev. D 16 (1977) 1791 
S. Weinberg, \Phys. Rev. Lett. 40 (1978) 2gf ; F. Wilczek, \Phys. Rev. Lett. 40 (1978) 279 



Nucl. Phys. B 166 (1980) 493 



J. E. Kim, Phys. Rev. Lett. 43 (1979) 103; M. A. Shifman, A. Vainstein and V. I. Zakharov, 



M. Dine, W. Fischler and M. Srednicki, [Phys. Lett. B 104 (1981) 199t A. P. Zhitnitskii, Sov. 
J. Nucl. 31 (1980) 260. 



L. F. Abbott and P. Sikivie, Phys. Lett. B 120 (1983) 133; J. Preskill, M. Wise and F 



Wilczek, Phys. Lett. B 120 (1983) 127; M. Dine and W. Fischler, Phys. Lett. B 120 (1983) 



137; M. Turner, Phys. Rev. D 33 (1986) 889; K. J. Bae, J. H. Huh and J. E. Kim 



JCAP0809 (2008) 005; L. Visinelli and P. Gondolo, Phys. Rev. D 80 (2009) 035024 



H. P. Nilles and S. Raby, Nucl. Phys. B 198 (1982) 102 



For recent reviews of axion/axino dark matter, see F. Steffen, Eur. Phys. J. C 59 (2009) 557 
L. Covi and J. E. Kim, New J. Phys. 11 (2009) 105003. 



T. Got o and M. Yamaguchi, [Phys. Lett. B 276 (1992) 103[ E. J. Chun, J. E. Kim and H. P. 
Nilles, \Phys. Lett. B 287 (1992) 123 . 



K. Rajagopal, M. Turner and F. Wilczek, [Nucl. Phys. B 358 (1991) 447 



E. J. Chun and A. Lukas, phys. Lett. B 357 (1995) 43 



L. Covi, J. E. Kim and L. Ros zkowski, [Phys. Rev. Lett. 82 (1999) 418$ L. Covi, H. B. Kim, 
J. E. Kim and L. Roszkowski, J. High Energy Phys. 0105 (2001) 033 . L. Covi, L. Roszkowski 
and Small, \j. High Energy Phys. 0207 (2002) 023 . 



A. Freitas, F. Steffen, N. Tajuddin and D. Wyler, [Phys. Lett. B 679 (2009) 270 



L. Covi, L. Roszkowski, R. Ruiz de Austri and M. Small, [J. High Energy Phys. 0406 (2004) 



003; A. Brandenburg, L. Covi, K. Hamaguchi, L. Roszkowski and F. Steffen, Phys. Lett. B 



(2008) 016 



617 (2005) 99; K. Y. Choi, L. Roszkowski and R. Ruiz de Austri, J. High Energy Phys. 0804 



H. Baer, A. Box and H. Summy, J. High Energy Phys. 0908 (2009) 080 



A. Chamseddine, R. Arnowitt and P. Nath, Phys. Rev. Lett. 49 (1982) 970| ; R. Barbieri, 
S. Ferrara and C. Savoy, phys. Lett. B 119 (1982) 343fc N. Ohta, Prog, fheor. Phys. 70 
(1983) 542; L. Hall, J. Lykken and S. Weinberg, phys. Rev. D 27 (1983) 235S. 



- 21 - 



[18] P. Sikivie and Q. Yang, \Phys. Rev. Lett. 103 (2009) 111301 



[19] H. Baer, S. Kraml, S. Sekmen and H. Summy, J. High Energy Phys. 0803 (2008) 056 ; H 



Baer and H. Summy, \Phys. Lett. B 666 (2008) 5| ; H. Baer, M. Haider, S. Kraml, S. Sekmen 
and H. Summy, JCAP0902 (2009) 002. 



[20] J. Ellis and K. Olive, \Phys. Lett. B 514 (2001) 114 . 

[21] F. Paige, S. Protopopescu, H. Baer and X. Tata, |hep-pti/0312045| ; 
http://www.nhn.ou.edu/~isajet/ 

[22] H. E. Haber, R. Hempfling and A. Hoang, \Z. Physik C 75 (1996) 539) . 



[23] D.Pierce, J. Bagger, K. Matchev and R. Zhang, \Nucl. Phys. B 491 (1997) 3 . 
[24] H. Baer, C. Balazs, A. Belyaev, \J. High Energy Phys. 0203 (2002) 042] . 



[25] H. Baer, C. H. Chen, M. Drees, F. Paige and X. Tata, \Phys. Rev. Lett. 79 (1997) 986 



[26] K. L. Chan, U. Chattopadhyay and P. Nath, Phys. Rev. D 58 (1998) 096004 ; J. Feng, 



K. Matchev and T. Moroi, Phys. Rev. Lett. 84 (2000) 2322 and Phys. Rev. D 61 (2000) 



075005] see also H. Baer, C. H. Ch en, F. Paige and X. Tata, \Phys. Rev. D 52 (1995) 2746j 



and \Phys. Rev. D 53 (1996) 624l| ; H. Baer, C. H. Chen, M. Drees, F. Paige and X. Tata, 



Phys. Rev. D 59 (1999) 055014; for a model- independent approach, see H. Baer, 



T. Krupovnickas, S. Profumo and P. Ullio, J. High Energy Phys. 0510 (2005) 02C. 



[27] J. Ellis, T. Falk and K. Olive, Phys. Lett. B 444 (1998) 367; J. Ellis, T. Falk, K. Olive and 



M. Srednicki, Astropart. Phys. 13 (2000) 181; M.E. Gomez, G. Lazarides and C. Pallis, Phys 



Rev. D 61 (2000) 123512j and [Phys. Lett. B 487 (2000) 313; ; A. Lahanas, D. V. Nanop oulos 



and V. Spanos, \Phys. Rev. D 62 (2000) 023515| ; R. Arnowitt, B. Dutta and Y. Santoso, [Nucl 



Phys. B 606 (2001) 59; see also Ref. 



[28] R. Arnowitt and P. Nath, Phys. Rev. Lett. 70 (1993) 3696; H. Baer and M. Brhlik, Phys 



Rev. D 53 (1996) 597; A. Djouadi, M. Drees and J. Kneur, Phys. Lett. B 624 (2005) 60 



[29] H. Goldberg, \Phys. Rev. Lett. 50 (1983) 1419 ; J. Ellis et al. \Nucl. Phys. B 238 (1984) 453 
P. Nath and R. Arnowitt, ^hys. Rev. Lett. 70 (1993) 3696|; H. Baer and M. Brhlik, Ref. 



V. Barger and C. Kao, Phys. Rev. D 57 (1998) 3131 



[30] M. Drees and M. Nojiri, [Phys. Rev. D 47 (1993) 376| ; H. Baer and M. Brhlik, \Phys. Rev. D 



57 (1998) 567| ; H. Baer, M. Brhlik, M. Diaz, J. Ferrandis, P. Mercadante, P. Quintana and 



Phys. Rev. D 63 (2001) 015007; J. Ellis, T. Falk, G. Ganis, K. Olive and 



X. Tata, 

M. Srednicki, \Phys. Lett. B 510 (2001) 236|; L. Roszkowski, R. Ruiz de Austri and T. Nihei, 



J. High Energy Phys. 0108 (2001) 024; A. Djouadi, M. Drees and J. L. Kneur, J. High 



Energy Phys. 0108 (2001) 055| A. Lahanas and V. Spanos, \Eur. Phys. J. C 23 (2002) 184 



[31] J. E. Kim and H. P. Nilles, \Phys. Lett. B 138 (1984) 150| . 



[32] K. Kohri, T. Moroi and A. Yotsuyanagi, \Phys. Rev. D 73 (2006) 123511] ; for an update, see 

see 



M. Kawasaki, K. Kohri, T. Moroi and A. Yotsuyanagi, Phys. Rev. D 78 (2008) 065011 



also J. Pradler and F. D. Steffen, \Phys. Lett. B 648 (2007) 224 . 
[33] S. Soni and H. A. Weldon, \Phys. Lett. B 126 (1983) 215| ; V.Kaplunovsky and J. Louis, \Phys 



Lett. B 306 (1993) 269; A. Brignole, L. Ibanez and C. Munoz, Nucl. Phys. B 422 (1994) 125. 



[34] For a recent analysis, see M. Carena, G. Nardini, M. Quiros and C. Wagner, Nucl. Phys. B 



812 (2009) 243. 



- 22 - 



[35] A. A. Affolder et al. [CDF Collaboration], Phys. Rev. Lett. 84 (2000) 5273; 



M. Fukugita and T. Yanagida, Phys. Lett. B 174 (1986) 45; M. Luty, Phys. Rev. D 45 



(1992) 455| ; W. Buchmiiller and M. Plumacher, \Phys. Lett. B 389 (1996) 73| and |/nt J. Mod 



Phys. A 15 (2000) 5047; R. Barbieri, P. Creminelli, A. Strumia and N. Tetradis, Nucl. Phys 



B 575 (2000) 61; G. F. Giudice, A. Notari, M. Raidal, A. Riotto and A. Strumia, Nucl 



Phys. B 685 (2004) 89; for a recent review, see W. Buchmiiller, R. Peccei and T. Yanagida, 



Ann. Rev. Nucl. Part. Set. 55 (2005) 311 



[37] W. Buchmuller, P. Di Bari and M. Plumacher, Annal. Phys. 315 (2005) 305. 



G. Lazari des and Q. Shan, \Phys. Lett. B 258 (1991) 305| ; K. Kumekawa, T. Moroi and T. 
Yanagida, \Prog. Theor. Phys. 92 (1994) 431 ; T. Asaka, K. Hamaguchi, M. Kawasaki and T. 
Yanagida, \Phys. Lett. B 464 (1999) 12|. 



[39] M. Ibe, T. Moroi and T. Yanagida, \Phys. Lett. B 620 (2005)1 . 
[40] I. Affleck and M. Dine, Nucl. Phys. B 249 (1985) 361. 



[41] H. Murayama and T. Yanagida, Phys. Lett. B 322 (1994) 349; M. Dine, L. Randall and S. 



Thomas, Nucl. Phys. B 458 (1996) 291 



[42] K. Jedamzik, M. LeMoine and G. Moultaka, JCAP0607 (2006) 010. 

[43] A. Brandenburg and F. StefTen, JCAP0408 (2004) 008. 

[44] R. Kolb and M. Turner, The Early Universe (Addison- Wesley, 1990). 

[45] For a recent review of direct, indirect and collider detection of neutralino dark matter, see H. 
Baer, E. K. Park and X. Tata, New J. Phys. 11 (2009) 105024. 



L. Duffy et al, Phys. Rev. Lett. 95 (2005) 091304 and Phys. Rev. D 74 (2006) 012006 



for a 



review, see S. Asztalos, L. Rosenberg, K. van Bibber, P. Sikivie and K. Zioutas, Ann. Rev 
Nucl. Part. Set. 56 (2006) 293j . 



[47] For some discussion of electroweak fine-tuning in supersymmetric models, see e.g. J. Ellis, K. 
Enqvist, D. V. Nanopoulos and F. Zwirner, Mod. Phys. Lett. A 1 (1986) 57; R. Barbieri and 



G. F. Giudice, Nucl. Phys. B 306 (1988) 63 



G. Anderson and D. Castano, Phys. Lett. B 



347 (1995) 300| and \Phys. Rev. D 52 (1995) 1693|; J. L. Feng, K. Ma tchev and T. Moroi 



Phys. Rev. Lett. 84 (2000) 2322| and \Phys. Rev. D 61 (2000) 075005 



[48] N. Arkani-Hamed, A. Delgado and G. F. Giudice, Nucl. Phys. B 741 (2006) 108; H. Baer, 



A. Mustafayev, E. K. Park and X. Tata, JCAP 0701 (2007) 017 and \J. High Energy Phys 



0805 (2008) 058 



- 23 - 



